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Abstract 

A series of calculations of a jet in crossflow 
usinq a three-dimensional finite difference model 
of the Navier-Stokes equations were made on a series 
of progressively finer grids. With a reasonable 
number of grid points (40x30x20), calculated jet 
penetration and mixing characteristics were found to 
compare favorably with the experimental measurements 
of Crabb, Durao, and Whitelaw. The calculated re- 
sults were grid dependant (i.e., numerical or false 
diffusion was present In the result). A finer grid 
calculation indicated thut turbulence model defi- 
ciencies may become more noticeable in the calcula- 
ten results as the magnitude of numerical diffusion 
is reduced. 

Introduction 

With the increasing costs of combustor develop- 
ment testing, a great deal of interest has focused 
on the use of numerical models to screen design 
changes or develop new combustor concepts. This 
t.vpe of corrputer-based design process now appears 
to be feasible through the use of three-dimensional 
combiistor performance models. Ultimately* 
these models may be used to greatly improve the 
durability and reliability of gas turbine 
combustors. 

A number of major restrictions must be overcome 
before this type of design methodology can be adop- 
ted. First, the proper physics must be incorporated 
into the differential equations usei’ in the combus- 
tor model. Second, numerical methods must accu- 
rately solve these differential equations. Finally* 
the accuracy of the resulting code must be assessed 
against fundamental data* and improvements must be 
made to the code to alleviate identified 
deficiencies. 

Currently available three-dimensional combustor 
performance models have yet to be thoroughly as- 
sessed. A few comparisons. have been made with 
actual combustor hardware, but these have not 
been conclusive. A more logical first step is to 
examine the extent to which three-dimensional hydro- 
dynamic processes can be calculated. One flow field 
of this type for which a great deal of experimental 
data exist is jets in crossflow. 

Jets in crossflow are particularly relevant to 
the gas turbine combustor designer. Cooling air 
jets (dilution jets) are used to control the hot- 
gas temperature profile entering the turbine. They 
are also used to set up aerodynamic patterns within 
the combustor, (which promote mixing and control lo- 
cal burning zone stoichiometry). As a result, the 
jet penetration and mixing characteristics of jets 
In crossflow are of primary concern in the combustor 
design process. 

There have been a number of previous calcula- 
tions of jets In crossf low,(^«^*®); however, 
these studies have been limited by two main factors. 


First, although a great deal of experimental data 
exist, rarely have important parameters such as the 
turbulence field, inlet velocity profiles, and jet 
mixing characteristics been fully measured. This 
limits the flow field quantities one can conpare 
and imposes the need to assume inlet boundary con- 
ditions for the calculation. Recent measurements 
greatly reduce this problem. Second, core 
storage and economy requirements have limited pre- 
vious calculations to coarse grid systems. This 
results in some numerical error being present in 
the computed solution, which can possibly call into 
question any conclusions drawn from these studies. 

The present report expands on this previous 
work by employing a series of progressively finer 
grid systems to calculate the single jet in cross- 
flow experimentally measured by Crabb, Durao and 
Whitelaw. These experimental measurements pro- 
vide a fairly complete collection of velocities* 
turbulence intensities, and jet concentration pro- 
files with measurements of the inlet field. The use 
of a series of progessively finer grid systems al- 
lows a differentiation between numerical errors and 
the hydrodynamic modeling assumptions embodied in 
the three-dimensional combustor code.^^ The re- 
sults of this comparison will provide additional In- 
sight into the deficiencies of the turbulence model 
and the code numerics. 

Mathematical Formulation 

Symbols 

U^ Mean velocity 

u^ Fluctuating velocity about mean 

U^ Mainstream velocity 

Re* Cell Reynolds number ■ U ax/u 

^ 2 2 ? 

K Turbulence kinetic energy » l/2(u + v ♦ w ) 

c Turbulence energy dissipation rate 

e/5. Normalized mean jet fluid concentratiun level 

*1 at jet orifice; ■ 0 in main stream be- 
fore jet 

Finite Difference Mode' 

The finite difference model employed in this 
study is In construction to those previously 

reported.' It provides the capability to 
analyze steady-state,' three-dimensional, elliptic, 
turbulent. Incompressible flow fields. Only the 
pertinent features of the CLde will be reviewed here 
as the’ literature'^»^) Is available for further 
details. 

The governing equations Include 


Continuity 
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Moment UiTi 
. 



Convection 


_i l£_ ♦ (2) 


Pressure 

gradient 


Diffusion 


Scalar Transport 


II iJL i Li- ii-il tt-iN 

3x^ * p ax^ “ 


Convection 


V- 

Diffusion 


( 3 ) 


Turbulence Model 



whore K and c are the kinetic energy and dissipa- 
tion rate of turbulence energy defined by 



whe 'e capitals indicate mean quantities and small 
let ers indi.:ite fluctuating quantities. The above 
pqu. tions correspond to the turbulence model of 
Jont s and Launder.^®) The constants used in the 
turbulence modol include 


« 0.09, - 1.44, C,2 - 1.92, 

0 |^ « 1.0, Oj. ■ 1.3, and * 0.9 


Solution Procedure 


Boundary Conditicns . Inflow conditions are of 
the Oirchlet type witiTthe main flow introduced 
approximately two and one-half jet diameters up- 
stream of the jet injection point. Jet boundary 
conditions were interpolated from the measure- 
nents(') unless otherwise stated. The rms tur- 
bulence intensity was chosen from experimental mea- 
surements as a uniform 7 percent of the jet velocity 
and the inlet length scale was based on 0.03 of the 
jet diameter. The downstream boundary condition was 
located approximately 14 jet diameters downstream of 
the jet injection point. 

Cross stream boundary conditions in the z di- 
rection (see Fig. 1) were modified from the original 
code to allow synmetry at the end points. The ori- 
ginal code provided for cyclic boundary conditions 
at the end points, which would have required the use 
of twice as many grid points to analyze a jet in 
crossflow. The provision of symmetry permitted the 
calculation to be bounded by the jet centerline. 

Grid Refinement . Hybrid differencing of the 
convective terms Tsee eq. (2) in the governing equa- 
tions requires that some form of grid refinement be 
used to obtain reasonably accurate solutions. Hy- 
brid differencing is a scheme whereby the type of 
difference formula employed is dependent on the 
Reynolds number of the computational cell (Rec). 
Second-order accurate central differencing is used 
when |Rec|< 2, rirst-order accurate upwind dif- 
ferencing is used when [Re^. | > 2, For calculations 
of practical concern to corrbuster designers, it is 
impractical to add enough grid points to maintain a 
Rec limit of 2, As a consequence, the computational 
domain is predominantly upwind differenced resulting 
in a significant truncation error. The order of 
this truncation error is similar to the diffusive 
terms in the qoverning equations, hence, the name 
"numericar or "false diffusion," To minimize the 
contribution of this numerical error term to the 
computed solution, the grid was refined employing 
the method previously demonstrated and briefly sum- 
marized here.'^^J 


Taking a Taylor series expansion of each con- 
vective term in the axial momentum equation, the 
truncation series can be^evaluated. For example, 
the convective term aU^/ax has a truncation 
error of 


T 


1 


<x + 


1 dV 

iT j 6 

Ja.2n 


The qoverning equations are represented by fin- 
ite difference approximations on a staggered grid 
system. Hybrid differencing of the convective 
terms was used to obtain stable, wiggle-free solu- 
tions, Central differencing was used for the other 
terms in the modelled equations. Each iteration of 
the solution begins with the finite difference forms 
of the momentum equation being sequentially solved. 

A Poisson pressure equation is then used to compute 
the pressure field. This in turn is used to "cor- 
rect" the velocity field to maintain continuity. 
Following this, the equations for K, c and • (in 
this case equals normal, zed jet fluid concentration) 
are solved using the corrected velocity field. Iter- 
ation on the computed flow field Is continued until 
the sum of the mass residuals is a low value, typi- 
cally less than 1 percent. This comprises the 
SIMPLE (Semi Implicit Pressure Linked Equations) 
algorithm as developed by Spaulding, 


where n Is an integer. If the initial term in tils 
series is used as a measure of the numerical diffu- 
sion, the grid can be adjusted to minimize the 1r- 
fluence of this false diffusion term in the solution 
of the overall equation. This* is done by making a 
preliminary calculation of the flow field, and from 
this computed field the Initial term In the trunca- 
tion series is approximated using central differen- 
cing. The net truncation error Is then obtained bv 
adding together the initial truncation error for all 
convective terms. The absolute value of this net 
convection truncation error is then compared with 
the magnitude of the other terms (pressure gradient 
and net convection) in the equation, which are also 
evaluated using central differences. This corrpari- 
son should then be made at each grid point in the 
three-dimensional domain. In theory, additional 
grid points are then added until the net contribu- 
tion of the truncation terms is small compared with 
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the other terms In the equation. This would then 
result In a grlrt-lndependant solution, given that 
the higher order truncation terms tend toward zero 
and that this procedure Is applied to each equation 
solved. More precisely, the solution will be 
spcond-ordcr accurate. 

This alone does not guarantee an accurate finite 
difference solution of the governing equations. 
Castrol^^I has shown that when both upwind and 
centi'al differencing are applied to the solution of 
• the one-dimonslonal convectlorv-dlffuslon equation, 
at high cell Reynolds numbers (|Recj > 2) upwind 
differencing Is more accurate. In fact. It Is this 
rationale has led to the use of hybrid dlf- 
> ferencing, This can be highly misleading, 

however, when eKtended to two or thre^dlmenslonal 
flows.^^^) The only Instance where hybrid dif- 
ferencing may be valid Is when the grid structure 
can be alined directly along flow streamlines. Ad 
hoc attempts to aline the rectangular grid system 
to flow streamlines are not likely to result In 
Improved computational accuracy. 

Aside from the question of whether or not 
second-order accuracy is better than first-order 
accuracy. It Is apparent that if the real effects 
of the turbulence closure are to be observed In the 
computed solution, the contribution of numerical 
diffusion mist be minimized, Cthe^'wise the physical 
diffusion 0^ the turbulence model may be overwhelmed 
by numerical diffusion contributed by the truncation 
error. For this reason, the grid refinement proce- 
dure was enployed. Unfortunately, enough grid 
points could not be added to the three-dimensional 
flow field to demonstrate a second-order accurate 
solution; however, (as will be noted In the section 
Results) numerical diffusion could he reduced enough 
so that the physical deficiencies of the turbulence 
model could be partially discerned. 

The use of a rational qnd refinement scheme is 
espec'^lly important In three-dimensional flows 
where the luxury of doubling the number of grid 
points in all directions cannot be afforded both 
from the standpoint of excessive computing times 
and core storage requirements. The grid refinement 
procedure has the advantage of demonstrating In 
which areas of the flow field the grid can most 
effectively be refined to minimize false diffusion. 

Flow Configuration , A schematic drawing of the 
free jet coordinate system Is given in Fig. 1, The 
X (axial) direction Is oriented In line with the 
mainstream flow. The free jet Is Injected along 
the y direction and normal to the main-stream flow. 
The cross stream direction Is along the z axis. 

Calculations were made for the flow conditions 
corresponding to the experimental measurements of 
. Crabb, Durao', and Whitelaw.v') The main-stream 

flow was specified at a uniform 12 m/s. The dilu- 
tion jet to main-stream velocity ratio was 2.3. 
Boundary conditions In the z and y direction were 
I located far enough away from the Jet orifice to 

avoid having any effect on the computed flow field. 

Results and Discussion 

In the following discussion, experimental versus 
calculated results are shown for three different 
grid systems: A coarse grid (20x20x12 grid points), 
a medium grid (40x30x20), and a fine grid (00x40 
x22). To determine where the grid should be r^ 


fined, the grid refinement procedure mentioned 
earlier was used. The purpoie of this examination 
of various finite difference grid systems Is two- 
fold, First, the use of several grid systems In- 
sures that conclusions drawn from the calculated 
results are not anomalies of the number of grid 
points used. Second, it Is Important to quantify 
the trade-off involved between numerical accuracy 
and computational expense. Absolute precison may 
not be warranted In englneerlnq calculations, where 
the Inlet parameters to the calculation may be un- 
certain and only some gross features of the flow 
f'ld may be of interest. But the converse danger 
of wSing too coarse a grid system must be avoided, 
otlcrwise Important flow field characteristics may 
be masked. By using several grid systems, the 
trade-off between accuracy and computational ex- 
pense is illustrated. 

The main factor to be observed In the three 
different grid calculations Is the extent to which 
numerical or false diffusion affects the calculated 
results. Other factors such as Inlet boundary con- 
ditions and convergence criteria were maintained 
constant for each calculation. 

Coarse Grid . Comparison between a coarse grid 
calcuTatTon of the free jet flow field and the ex- 
perimental measurements! is shown In Fig. 2. 

The grid system, 20x20x12 (number of points In the 
X X Y X 2 directions), is completely Inadequate to 
accommodate the sleep velocity gradients In the 
flow field. The axial velocity peak evident In the 
experimental measurements at Y/0 ■ 2.b Is essen- 
tially missing In the computed flow field (see Fig. 
2(a)). The same effect is evident further down- 
stream in Fig. 2(b), but the general S-shaped 
nature of the flow is fairly well-predicted. The 
profiles of normalized jet fluid coincide fairly 
well but the calculate^ peak Indicates that the jet 
is slightly over-penetrating (see Fig. 2(c))- 

A gualitatlve examination of the cross-stream 
vortex formed In the course grid calculation Is 
displayed by velocity vector plots (see Fig. 3), 

The vortex structure Is just being formed at a 
cross section bisecting the jet orifice (see 
Fig, 3(a)), Just downstream of the jet orifice, 
the vortex can be seen to have gained strength (see 
Fig, 3(b)), The center of the vortex has also moved 
slightly closer to the Jet centerline (at 2/D • 0) 
The vortex can be seen to have expanded In size 
while the center Is slightly displaced outward (see 
Fig. 3(c)). In a qualitative sense, this vortex 
structure agrees with the measurements of Ref, 7. 

The loci of the vortex center could be traced from 
the experimental measurements to be drawn back Into 
the wake region behind the jet, and then be trans- 
ported outward further downstream. However, the 
precise location of this vortex center and Its 
movement Is not well-predicted. For example, the 
experimental measurements Indicate that the outward 
movement should occur around X/D « 3,25. The calcu- 
lations clearly predict this will occur much closer 
to the jet orifice, 

A consequence of this strong, calculated vortex 
structure can be seen In Fig. 4, Here the experi- 
mental and calculated normalized Jet concentration 
profiles at an axial location of X/D • 6 are com- 
pared, The strong cross stream vortex causes the 
ft:ax1mum of the Jet concentration to occur off the 
Jet centerline. This physically agrees with the 
experimental results; however, the strong cross 
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Stream vortex may also augment the penetration of 
jet fluid in the Y direction with the net result 
leaolng to a poor comparison between calculated and 
experimental results. 

Medium Grid . Aaditional grid points were added 
to the computational domain to observe whether any 
improvement over the coarse grid calculation could 
he obtained. A comparison between the calculated 
and experimental results can be seen In Fig, 5. 

The use of 40x30x20 grid points provided a greatly 
improved accommodation of the velocity gradients In 
the flow field. Axial velocity profiles at two 
different axial locations display an Improved com- 
parison between calculated and experimental results 
{Fiqs. 5(a) and (b)). In general, the jet penetri>- 
lion appears to coincide well with the experimental 
measurements. In Fig. 5(c) the axial turbulence 

intensity (^u^) is compared at an axial location 
of X/D « 2, The magnitude of the axial turbulence 
intensity compares favorably with the experiment 
around Y/0 « 1.5, but is underpredicted elsewhere. 

The moderately good agreement shown in Fig 5(c) 
can be misleading. Explicit data are not available 
at this axial location, but the measurements of 
Ref. 7 indicate that in the jet wake region, the 
magnitude of the cross stream (2 direction) turbu- 
lent fluctuations is much greater than the magni- 
tude of the axial fluctuations (X direction). The 
i?otropic turbulence model used in the calculation 
procedure assumes an equal magnitude of turbulence 
ii'tensUy in all directions; therefore, the cross- 
stream fluctuations are likely to be significantly 
underpredicted by this model. This is a common 
problem with the isotropic assumption and has been 
noted elsewhere. 

In Fiq. 5(d) profiles of axial turbulence 
intensity and normalized jet fluid concentration 
are shown at an X/D of 6. Here, the turbulence 
intensity is consistently underpredicted and the 
sane isotropic turbulence limitation is likely to 
bt true. However, the isotropic assumption should 
br less restrictive since the flow will have a ten- 
d(ncv toward isotropy far downstream of the jet 
01 if ice. 

One of the interesting experimental findings of 
Ref. 7 is that the profiles of mean axial velocity 
and jet fluid concentration uo not coincide. The 
frt e-stream fluid is accelerated around the jet and 
cai ses the velocity peak around Y/0 - 3.5 (see Fig, 
5(t). The concentration of jet fluid peaks around 
Y/f * 2.8 (i.e., in the wake region). Expei imen- 
ta ly, the jet fluid concentration profiles were 
fO' nd to coincide more appropriately with the axial 
tui tulence Intensity profile, instead of the axial 
velocity profile. Both trends are well-reproduced 
in the calculated results (Figs. 5(b) and (c)). 

Further verification of the accurate prediction 
of jet penetration can be seen in Fig. 6, Calcu- 
lated velocity vectors are displayed against an 
enpirical correlation for jet penetration, 

The velocity vectors match the empirical trajectory 
quite well. Another positive feature of the com- 
puted flow field can be seen In this figure. The 
upwelling, or positive V velocity, that can be seen 
In the recirculation zone behind the jet qualita- 
tively agrees with the results of Ref. 7. 


A less favorable aspect of the conputed flow 
field can be seen In the velocity vector plot of 
Fig. 7. Here the y-z plane at an axial location of 
x/0 • 2 displays an unusual vortex structure. Two 
coflowing vortices are embedded In the overall 
cross-stream vortex. This structure, which was not 
apparent In the coarse grid results (Fig. 3), cannot 
be supported by the experimental measurements. It 
appears i,o be a nonphysical artifact of the calcu- 
lation, which may be related to the Isotropic tur- 
bulence model. If the calculated cross stream tur- 
bulence intensity were greater, these two vortices 
might collapse into one large cross-stream vortex. 

Despite this nonphysical aspect of the velocity 
field, an excellent agreement exists between calcu- 
lated and experimental results for the normalized 
jet fluid concentration profiles at X/D - 8 (see 
Fig. 8). The most surprising feature of this agree- 
ment Is the accuracy of the cross stream (z direc- 
tion) transport of jet fluid. As noted earlier, 
the cross stream fluctuations could be Inferred to 
be larger in magnitude than the axial turbulent 
fluctuations in the jet wake region. Therefore, 
the inadequacy of the isotropic turbulence mode^ 
should have exhibited itself as a lack of turbulent 
diffusion In the z direction. The reason for this 
discrepancy is a result of numerical or false dif- 
fusion, as will be made clear In the next series of 
calculations. 

results of a fine grid calcula- 
tion (95x30x?2 ) are compared against experimental 
measurements in (^ig, 9). Here the addition of 
grid points has slightly improved the axial velocity 
comparisons (Fiq. 9(a)) while lessening the agree- 
ment on jet concentration profiles (Fig. 9(b)). 

The recirculation zone velocity is matched more 
closely than the previous results and the velocity 
peak at Y/D ■ 2.4 is sliqhtly closer to the experi- 
mental value (see Fiq. 9(a)), The calculated jet 
concentration profiles indicate a lack of jet 
spreading in the Y direction. This is especially 
noticeable in the region around Y/D - 3 - 4, The 
0,20 contour does not spread in the Z direction as 
far as previous results. The maximum penetration 
is about Z/D ■ 1.0 as opposed to an experimental 
value of about 1.3 (1,3 is also the maximum pene- 
tration determined In the medium grid calculation). 
The O.lO contour does spread correctly in th» 2 
direction. 

The results (Fig, 9) are still grid dependent 
and the reason for this is seen In Fig. 10, Three 
of the terms in the axial momentum equation 
(eq. (2)) are shown at an axial location of X/D 
• 1,3 near the jet centerline. (The diffusion term 
is excluded for clarity). These terms were eva- 
luated using central differencing of the computed 
flow field In a manner consistent with the proce- 
dure defined In the section Grid Refinement. The 
main feature of the Indicated profiles is that the 
convection truncation error (numerical or false 
diffusion term) Is negligible for much of the flow 
field but becomes quite significant around a Y/0 
of 2, The axial velocity profile (Fig. 10(b)) 
corresponds to the same location as that in 
Fig, 10(a), The truncation errors can be ieen to 
correspond to the velority peak of the fluid acce- 
lerating around the jet. This Is a red on of very 
steep velocity gradients, which causos’the locally 
large truncation errors. The dtvt>lopment of the 
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velocity pe«^k is largely 
diffusion {f ios. 10(a) and (b)). 

Truncation errors are even more serious In the 
solution of the scalar transport equation (eq. 3). 
(This equation is solved to determine jet fluid con- 
centration profiles - Fig, 9(b). In the solution 
of this equation, the convection terms (determined 
by hybrid differencing) are set equal to the diffu- 
sion terms (determined by central differencing). 

There is no pressure gradient term. This results 
in the convection truncation error having a large 
impact on the solution of the equation. For exam- 
ple, Fiq. 11 displays the magnitude of the convec- 
tion and truncation terms in the scalar transport 
equation at X/D » 1.9 and Z/0 « 0.9. The trunca- 
tion terms contribute significantly to the computed 
result. This explains how the jet fluid concen- 
tration profiles can spread in a manner not con- 
sistent with the physics of the turbulence model. 

The solution is so badly swamped by numerical dif- 
fusion in some regions of the flow field that it 
almost controls the computed results. It is only 
in the fine grid results (Fig. 9(b)) that an indi- 
cation is given of what an error-free solution might 
vield. The concentration profiles would not spread 
quite so far in the Y or Z direction as was shown 
in the 0.20 concentration profile. The agreement 
seen in the 0.10 profile is simply a result of the 
lack of grid points in this region of the flow field. 

Compu ti ng T imes . The coarse grid calculations 
displayed in this report required approximately 
iiO CPU min of computing time on an IBM model 
370-3033. The medium grid calculations required 
about 2 CPU hours and the fine grid calculations 
about 10 CPU hr. 

Summary . Comparisons between experimental and 
calciilaTod results for three different grid systems 
have been presented. All grid systems were subject 
to some numerical error (numerical diffusion). The 
coarse-qrid calculations exhibited a high degree of 
numerical diffusion resulting in a obscuration of 
the jet velocity peak, especially near the jet 
orifice. The medium-grid calculations exhibited 
noticeably less numerical diffusion with an unoer- 
prediction rf the magnitude of the jet peak velocity 
but a correct penetration trajectory. The jet con- 
centration profiles for this grid compared very 
well with experiment. The fine-grid calculation 
exhiuited a slightly better comparison of the jet 
peak velocities' but a worse correspondence between 
experimental and calculated jet concentration 
profiles. 

All calculations provided a qualitative predic- 
tion of some features of the single jet flow field. 
For the cross-stream vortex, the agreement became 
less qualitative as additional grid points were 
added, which is in physical agreement with the tur- 

• bulence model. 

Previous calculations of jets in cross flov/^) 
have indicated an excellent correspondence with 

• experimentally determined jet trajectories. It Is 

postulated that the pressure field essentially con- 
trols jet penetration and, hence, any inaccuracies 
in the turbulence closure. or the differencing scheme 
are of lesser importance.''' From the results of 
this report, it appears that this may or may not be 
true, depending on the number of grid points used 

in the calculation. Coarse grid calculations in- 
dicated some overpenetration of the jet, but with a 


was well predicted. 

Equally important are the jet fluid mixing 
characteristics; and these were found to be very 
sensitive to numerical and turbulence model Inac- 
curacies. Unfortunately, the response of these 
variables to grid refinement is only sluggish at 
best. 

Concluding Remarks 

What is illustrated quite clearly is the need 
for an improved numerical scheme to 
merical diffusion. A variety of techniques'^®*^'' 
may provide a resolution to this problem. However^ 
this is not to say that the current three-dimen- 
sional combustor models cannot be used to make 
engineering calculations. For the flow field in 
this investigation, the medium grid (40x30x20 grid 
points) provided a good compromise between computa- 
tional expense and accuracy. For other types of 
flows this would doubtless vary and currently can 
only be determined by a post analysis of the trun- 
cation error. 

Conclusions 

Employing a three-dimensional finite-difference 
model to analyze the jet flow field,'') the fol- 
lowing conclusions were determined: 

1. With a reasonable number of grid points 

( approximately 40x30x20), calculated jet pene- 
tration and concentration profiles agreed well 
with experimental measurements. 

2. For the cross- stream vortex, the agreement 
between experimental and calculated results 
became less qualitative as additional grid 
points were added, indicating a defiency of the 
isotropic turbulence model. 

3. The calculated results of the finest grid 
examined (90x40x22) were grid-dependent. An 
iirproved numerical scheme is required to remove 
the effects of numerical diffusion for the flew 
geometry examined. 
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(b) X/D-6. 
lc)X/0-6. 

Figure?. * Coarse grid (20x20x12) 
caicuiaticn of a single, free Jet 
Velocity profiles at a) XfD-2, b) 
X/D>^ and c) normalized jd 
fluid concentration profiles at 
XfD^ All profiles shoem are 
throu^ the jet centerline. 
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Figure 4. - Coarse grid 
(20 x2Ck 12) calculation of 
Jet fluid concentration 
profiles and comparison 
with experiment at X/D^. 
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(a)X/D-a 12$. 
IWX/D-a«5. 
(c)X/D-L375. 

Figure 3. - Velocity vector plots 
of a coarse grid GQk20k 12) cal- 
culation of the single Jet flow 
field displaying the cross- 
stream vortex. 
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(a) Axial valadty profiles along Jet centerline 
atX/D-2. 


(b) Axial velocity profiles along Jet centerline 
atXIO-6. 

(c) Turbulence profiles along Jet centerline 
atX/0-2. 

(dl Turbulence and Jet fluid concentration 
profiles along Jet centerline at X/D*b. 

Figure 5. - grid (4(k30K20l cal* 
culation of sik< 'e, free Jet flow field and 
comparison with «(periment 
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Figures. - Medium grid 
<409(30(20) alculaticn of 
let fluid concentration 
profiles and comparison 
with experiment atX/D«6. 



ia) Axial v«i(^ity profl«a at Xf 0-2. 

(b) Jot fluid contratlon profile at 
X/D-8. 


Figure 9. - Fine grid (90(40(22) cal- 
culation of a single )et 



ORIGINAL PA^c \'J 
Of POOR QUALITY 



(al Convection, pressure gradient and convec- 
tive truncation error at X/D-L 3 nea r the Jet 
centerline. 


(b) Axial velocity at X/D-L 3 near the Jet center- 
line. (Grid 90x40x22.) 

Figure 10. - Comparison of the terms in the 
axial momentum equation. 



Figure 11. - Comparison of convective 
^rms in the scalar transport equation 
at X/D-L 9 and Z/O-a « (Grid 90x40x22). 



